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We study a two-species bosonic Hubbard model on a two-dimensional square lattice by means of 
quantum Monte Carlo simulations. In addition to the usual contact repulsive interactions between 
the particles, the Hamiltonian has an interconversion term which allows the transformation of two 
particles from one species to the other. The phases are characterized by their solid or superfluid 
properties and by their polarization, i. e. the difference in the populations. When inter-species inter- 
actions are smaller than the intra-species ones, the system is unpolarized, whereas in the opposite 
case the system is unpolarized in even Mott insulator lobes and polarized in odd Mott lobes and also 
in the superfluid phase. We show that in the latter case the transition between the Mott insulator of 
total density two and the superfluid can be either of second or first order depending on the relative 
values of the interactions, whereas the transitions are continuous in all other cases. 
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I. INTRODUCTION 

Since their first experimental realization, atomic Bosc- 
Einstcin condensates have been used as a tool to study 
strongly interacting quantum systems. This made possi- 
ble the study of exotic phases, especially in systems gov- 
erned by Hamiltonians not easily realized in condensed 
matter. Soon after the initial work on spin-0 bosons 
much effort was directed at the problem of mixtures of 
different kind of particles: fermions [H or mixtures of 
bosons and fermions However, in these studies, the 
effective spin degrees of freedom are frozen and the num- 
ber of particles of each type is generally fixed. Purely op- 
tical trapping techniques allow to trap the atoms without 
freezing the internal degrees of freedom. These systems 
exhibit both magnetic and superfluid properties offering 
the opportunity to study the interplay between these two 
effects. Unlike in simple mixtures, spin-spin interactions 
are present in these systems and can be either ferromag- 
netic ( 87 Rb) or antiferromagnetic ( 23 Na)[J, 

Whereas current spin- 1 experiments are in continuous 
space, a simple model was proposed @ for a system on an 
optical lattice with only two internal low energy states. 
The proposed system consider neutral polarizable atoms 
with three degenerate atomic ground states and three 
degenerate excited states characterized by the magnetic 
quantum number S z = 0, ±1. In addition to generating a 
periodic (optical lattice) potential, the counter propagat- 
ing lasers couple these internal ground and excited states 
by V and A transitions which leads to only two low energy 
cigenstates denoted respectively and A. Such particles 
with only two internal effective degrees of freedom are 
referred to as spin-1/2 bosons. As in the spin-I case, 
the interaction between these spins can be ferromagnetic 
or antiferromagnetic and the presence of an optical lat- 
tice allows the system to become strongly correlated with 
superfluid (SF) or Mott insulator (MI) phases. In pre- 
vious work Q, we extended the original mean-field the- 



ory (MFT) approach Q and studied the one-dimensional 
system with exact quantum Monte Carlo (QMC) simula- 
tions. Related spin-1 models were considered using MFT 
@, [| or in one dimension using QMC 0. A similar 
model with two species of bosons that can convert into 
each other was also studied [lf| . 

If the spin-spin interactions for the original spin one 
bosons are ferromagnetic (F) , the resulting on-site repul- 
sions for spin 1/2 bosons are smaller between different 
particles than between identical particles [1, 0] ■ In that 
case MFT and QMC both show the system never polar- 
izes Q; in other words, the populations of the two species 
always remain equal. If the repulsion is strong enough, 
the system is in a Mott insulating phase for commensu- 
rate densities and is otherwise superfluid. The transitions 
between the superfluid and the Mott phases are always 
continuous. 

In the case of antiferromagnetic (AF) interactions for 
the original spin one bosons, the repulsion between dif- 
ferent particles is stronger than the repulsion between 
identical particles @, 0] in the resulting spin 1/2 model 
and the system exhibits a richer phase diagram. First, 
the Mott phases are polarized for odd densities while they 
are unpolarized for even ones. In addition, the superfluid 
phase is always polarized. Finally, MFT predicts second 
order transitions between the odd density Mott phases 
and the superfluid whereas the transitions between even 
MI and SF can be first or second order depending on the 
strength of the interaction. This was not observed in the 
one dimensional QMC simulations Q. This discrepancy 
was not surprising since MFT often fails in reduced di- 
mensionality and also since first order phase transitions 
are generally absent in one dimension. Another feature of 
the AF regime is a phase transition between a polarized 
and an unpolarized phase inside the p = 1 (p being the 
total particles density) Mott phase at low temperature 
observed with the QMC simulations. 

In this paper, we extend the QMC study of this sys- 
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tern to the two dimensional square lattice. The paper 
is organized as follows. In section II, we will review the 
model and the techniques used to study it. Sections III 
and IV will be devoted to the presentation of our results 
for the cases of larger repulsion between different or iden- 
tical particles, respectively. We summarize these results 
in Section V. 



II. SPIN-1/2 MODEL 

The model we will study has two species of bosonic 
particles governed by the Hamiltonian |6[ , 
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where a\ v creates a particle of "spin" a = 0, A on site 
r = [x, y) of a periodic square lattice of size L x L. The 
n av = <r aT a aT are the corresponding number operators 
counting the particles of type a on site r. The Hamil- 
tonian includes the conventional hopping term (Eq. ([T])) 
which plays the role of a kinetic energy term for lat- 
tice systems with the associated hopping parameter t 
that is used as the energy scale. To study the system 
in the grand canonical ensemble, a chemical potential 
term is added to the Hamiltonian ([1]). There are two 
density-density interaction terms (Eq. @). The first de- 
scribes repulsion between identical particles on the same 
site with an associated energy Uq > 0. The second de- 
scribes the on-site repulsion between particles of different 
types. Depending on the value of U2, this repulsion can 
be stronger (U2 > 0) or smaller (U2 < 0) than the repul- 
sion between identical particles. In this work, [ C/2 1 will 
remain smaller than Uq in order to keep only repulsive 
interactions. In most of our work we kept U2/U0 fixed as 
Uq is changed to map the phase diagram. Finally the 
last term of the Hamiltonian describes conversion be- 
tween the species: As two identical particles collide on 
the same site, they can be converted into two particles 
of the other kind. It was shown in Ref. @ that the ma- 
trix element for this conversion is equal to U 2 /2, that is, 
it is essentially the difference of interaction energies be- 
tween identical and different particles. The sign of the 
conversion term (Eq. fl3|_) can be chosen freely due to a 
symmetry of the model [7j . Here our convention has the 
opposite sign of the original paper Q. 

The only term in the Hamiltonian that couples dif- 
ferent sites is the hopping term in ([T]). MFT 0, @] iso- 
lates one site r and couples it to surrounding sites r' 
via the mean values of destruction or creation operators 



results in a one-site Hamiltonian which can be easily di- 
agonalized numerically. The ground state energy is then 
minimized with respect to the two mean-field order pa- 
rameters tpo and ipA ■ Normal or insulating phases are ob- 
tained when ipo = = whereas superfmid phases oc- 
cur whenever one of the tp is non zero. To study this sys- 
tem exactly, we used the stochastic Green function (SGF) 
quantum Monte Carlo algorithm [12, [l3j . This algorithm 
works in the canonical or grand canonical ensembles at 
finite inverse temperature /3 = 1/kT. We generally used 
/? = 2L/t which is usually a low enough temperature to 
obtain results that have converged to their ground state 
limit for a system of linear size L with moderate inter- 
actions. However (see below) for the strongly interacting 
regimes, it is sometimes necessary to use lower temper- 
atures (up to j8 = AL/t). In the canonical case with N 
particles, the chemical potential is defined as the discrete 
difference of the energy n(N) = E(N + 1) - E(N) which 
is valid in the ground state where the free energy is equal 
to the internal energy. The total particle density, p, is ei- 
ther fixed in the canonical case or fluctuates in the grand 
canonical one. Densities for particles of types and A are 
called po and p\, respectively, and are not conserved due 
to the conversion term Eq. §3§. The supcrfluid density is 
given by the fluctuations of the winding number [l4| 
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As explained in Ref. [15| . it is not meaningful to define 
individual supcrfluid densities for the and A particles 
as their numbers are not conserved separately due to the 
conversion term Eq. ([3|) . We also studied the single par- 
ticle Green functions, 
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where a — 0, A. The Fourier transforms of the Green 
functions, Eq.©, give the momentum distributions, 
Po-(k). The total density at zero momentum, p(k = 0), 
is given by 



P (k = 0) = 5> ff (k = 0). 
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Another useful quantity to characterize the Mott insu- 
lator is the two-particle anti-correlated Green function 
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which is a measure of the exchange process whereby, for 
example, a particle is annihilated at site r and a A 
particle at site r + R while a particle is created at site 
r + R and a A particle at site r. If perfect phase co- 
herence is established by means of particle exchange, G a 
reaches its limiting upper value of poPA at long distances 
R. 
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FIG. 1: (Color online) Zero temperature MFT phase diagram 
of the bosonic spin-1/2 model in the Ui > case. The dome- 
shaped surface shows the quantum phase transition between 
MI (below the surface) and SF regions (above). The grey 
areas indicate the regions where the transition is first order. 
Everywhere else the transitions are continuous. The transi- 
tions are also second order for V% = 0. 



III. POSITIVE U 2 CASE 



In the case where U2 > 0, MFT predicts at zero tem- 
perature [tJ that the system is either in a polarized su- 
per-fluid (SF) phase or in incompressible Mott insulat- 
ing (MI) phases. In the polarized superfluid phase, the 
symmetry between the two populations is broken and 
the density of one species becomes larger than the other. 
For total fillings which are odd multiples of the number 
of sites, i.e. p = 1, 3 ■ ■ • , the MI phases are also polarized 
whereas for even total densities, they are not: po = p\- 
The MI-SF transition is predicted to be continuous for 
odd density Mott phases. For even densities, the transi- 
tion is found to be first order at the tip of the MI lobes for 
small U2/U0 and becomes continuous as U2/U0 increases 
(see Fig. [T]). In the strongly interacting limit t/Uo 0, 
the Mott phases of odd density p are found for chemical 
potentials p in the interval p — 1 < p/Uo < p — U2/U0, 
whereas the even density Mott phases are observed for 
p — 1 — U2/U0 < p/Uo < p. Hence as observed in Fig. Q] 
or Fig. [21 the odd density Mott regions shrink as U2/U0 
is increased and disappear for U2/U0 = 1 in this t = 
limit. 
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FIG. 2: (Color online) QMC phase diagram at j3 = 2L/t 
or iL/t for U2/U0 = 0.1 and several linear system sizes L. 
The superfluid phase is always polarised and the p = 2 MI 
is not polarized due to the effect of the interspecies exchange 
term. For the p = 1 MI phase, the system is polarized for 
the larger values of t/Uo but then becomes unpolarized as 
t/Uo decreases at finite kT = t/2L. As T decreases further 
to kT = t/4L polarization appears (see Fig. |HJ. The dotted 
lines are the mean-field results. The region between the two 
MFT curves in the p = 2 lobe shows the coexistence zone due 
to the predicted first order transition. 



A. Phase diagram 



To map out the phase diagram, we employ the SGF 
algorithm in its canonical implementation. The chemical 
potential is calculated using energy differences to deter- 
mine the boundaries of the MI lobes for p = 1,2. In 
addition, we studied the histograms of the densities of 
the two species to determine whether a phase is polar- 
ized or not. We found results similar to those found in 
one dimension: The superfluid phase is always polarized, 
the p = 2 MI phase is not polarized and the p = 1 shows 
a transition between a polarized MI and an unpolarized 
MI as t/Uo decreases. This point will be discussed in 
more detail in Sec. IIII CI The resulting phase diagram 
is shown in Fig. [5] for several linear lattice sizes, L, and 
(3 = 1/kT = 2L/t or (3 = 4L/t. 



Using QMC simulations in the one-dimensional case 
, wc observed the polarized SF phase, the polarized 
MI lobes for p = 1 and the unpolarized MI for p = 2. 
However, contrary to MFT predictions, all transitions 
were continuous; as is often the case in one dimensional 
systems, no first order transitions were found. We also 
observed, at low temperature, that as t/Uo decreases (al- 
ways keeping U2/U0 constant) the p = 1 MI changes 
nature from polarized to unpolarized. 



The agreement between the QMC results and MFT 
is quite good. The Mott lobes obtained with MFT are 
shown as dotted lines in Fig. [5] and are close to the ex- 
act boundary for small values of t/Uo with disagreement 
increasing as the tips of the lobes are approached due to 
increased quantum fluctuations. As expected, the agree- 
ment is much better than in one dimension where we 
had found a factor of two difference between the MFT 
predictions and the observed lobe tips. 
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FIG. 3: The SF density, Eq. § as a function of t/U at 
fixed total particle density for U2/U0 — 0.1. The transition 
from the MI (small t/Uo) to the SF phase is signaled by p a 
becoming finite. For p = 1 (a) the transition is continuous 
whereas the jump in p a marks the presence of a first order 
phase transition for p — 2 (b). See Figs. [4] E] and[6] The inset 
shows the transition region and the jumps in more detail. 



FIG. 4: (Color online) The superfluid density from QMC and 
MFT calculations as a function of t/Uo for p = 2 and for 
several positive values of U2/U0. Both methods exhibit a 
jump in p s indicating the presence of a first order transition 
at the tip of the p = 2 Mott lobe for small but finite values 
of U2/U0. The jumps increases from at U2/U0 = to a 
maximum for U2/U0 — 0.05 and then decreases back to zero 
for U2/U0 > 0.25. 



B. Phase transitions 

An important difference between one and two dimen- 
sions is that first-order phase transitions may appear in 
the latter whereas they are essentially excluded in the 
former. Since MFT predicts a first order phase transi- 
tion for even Mott lobes, we used QMC to study this 
transition in great detail. 

Figure [3] shows the evolution of the superfluid density 
p s at the MI-SF transition at fixed total density while 
varying t/Uo, in other words the transition at the tip 
of the Mott lobe. We observe that in the p = 1 case 
(Fig. Eta)), the transition is continuous. There is no ev- 
idence at all for a first order transition for this case. On 
the contrary, in the p — 2 case (Fig.|3[b)), a jump in the 
superfluid density indicates a first order phase transition. 

Comparing this p = 2 case with MFT results for differ- 
ent values of U2/U0 in Fig.|U we observe very similar be- 
havior. The transition is continuous for U2/U0 = 0. Then 
it becomes discontinuous for small values of U2/U0 < 0.25 
and is again continuous for U%/Uq > 0.25. The jump in 
the superfluid density varies continuously from zero at 
U2/U0 = to a maximum for U2/U0 ~ 0.05 and then de- 
creases back to zero as U2/U0 is increased further. Figure 
[5] shows similar behavior for the particle density at zero 
momentum, p(k = 0), where, once again, the jump at 
the transition is observed for small values 0HJ2/U0. 

To confirm the presence of first order phase transitions 
near the tip of the MI lobe for even densities, we studied 
the behavior of p as a function of p, as one cuts across the 
lobe at fixed t/Uo using QMC simulation in the canonical 
and grand canonical ensembles. In the canonical enscm- 
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FIG. 5: (Color online) The density in the k = momentum 
p(k = 0) as a function of t/Uo for p = 2 and for several values 
of U2/U0. As for p s , Fig. [4] we observe a jump in p(k = 0) 
for small but finite values of U2/U0, confirming the presence 
of a first order transition. 



ble, a first order transition is signaled by negative com- 
pressibility [ll|, k cx dp/ dp, < 0. In the grand canon- 
ical ensemble, there will be a corresponding discontinu- 
ous jump in the p versus p curve. Figure [6] shows both 
these cases. In Fig.^a) the canonical simulations clearly 
show negative k just before and after the Mott plateau at 
p = 2. On the other hand, the grand canonical ensemble, 
Fig. IHJb), shows discontinuous jumps in p at the corre- 
sponding values of /i. The canonical and grand canonical 
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FIG. 6: (Color online) The density, p, as a function of the 
chemical potential, p, close to the tip of the p = 2 Mott lobe. 
The canonical simulation (a) shows a negative compressibility 
region, k oc dp/dp, < 0. The grand canonical simulation 
(b) exhibits a corresponding jump in the density. Thus both 
approaches show the presence of a first order transition. 



simulations are in quantitative agreement on the size of 
the unstable region which is extremely narrow; the sys- 
tem is stable for densities smaller than p = 1.95 or larger 
than p = 2.05 for the chosen value of t/Uo. 



C. Polarization and nature of the Mott phases 

We now analyse in more detail the polarization prop- 
erties of the Mott phases. In the p = 1 MI phase, 
Fig. [7] shows that a transition between unpolarized and 
polarized regions occurs for t/Uo ~ 0.02 for L = 8 and 
j3 = 2L/t. This is well inside the MI region as can be 
seen from the phase diagram, Fig. [5] This possibility 
of a polarization transition deep in the MI lobe requires 
closer examination. 

In the p — 1 MI phase, the system is frozen in a state 
with one particle per site. Fluctuations around this state 
will occur when particles hop around. However, events 
where two particles are converted to the other species 
are negligible since this requires double occupancies. Ne- 
glecting conversions and taking the hopping term as per- 
turbation to second order, the model can then be mapped 
onto an anisotropic Heiscnbcrg model 17 1 where the pres- 
ence of a A particle on a site corresponds to an up spin 
along the z axis and a particle corresponds to a down 
spin. The interactions between spins in the xy plane and 
along the z direction are given by [l7| 
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FIG. 7: (Color online) Histogram of the density of A particles 
for Ui > as a function of t/Uo- The total density is fixed 
at p = 1 and the histogram of the particles is the image of 
this distribution with respect to the line pA = 0.5. The other 
parameters are L = 8, /3 = 2L/t, and U2/U0 = 0.1. The 
system polarizes for t/Uo — 0.02, well inside the p — 1 Mott 
lobe. Thus at this low, albeit finite, temperature, there is a 
transition between a polarized Mott phase at large t/Uo and 
an unpolarized Mott phase at lower t/Uo. As shown in Fig. [8] 
the complete lobe is polarized at T = 0. 



For U2 > 0, we see that the couplings J xy and J z are 
always negative, i.e. the effective model spin interactions 
are always ferromagnetic. We also see that \J Z \ > \J X y\, 
which means that ferromagnetic order will develop along 
the z direction. In other words, the system will become 
polarized in terms of or A particles in the ground state 
limit. 

However, we also see that the energy associated with 
the polarization of the system is t 2 /Uo and becomes very 
small in the large Uo limit. On the other hand, the his- 
togram in Fig.[7]is at fixed temperature, (3 = 2L/t, which 
means that the system is no longer in the zero temper- 
ature limit for small values of t/Uo < 0.02 for L = 8. 
Then excitations of these spin degrees of freedom occur 
and the system is no longer polarized. 

This argument is confirmed by QMC simulations done 
at lower temperature: In Fig. [8] we show the evolution of 
the histogram of p\ for given values of Uo and U2 and for 
several inverse temperatures from /3 = 2L/t to /? = 4L/t. 
We observe that as the temperature is decreased, the 
system polarizes. We then conclude that, in the ground 
state limit, the entire p = 1 Mott phase is polarized, 
although one needs extremely low temperatures in order 
to observe the polarization for large Uo- 

Figure [9] shows single particle and anticorrelated Green 
functions in the MI (a) and SF (b) phases. As expected, 
we see in the Mott phase, Fig. E{a), that the individual 
Green functions Go and Ga decay exponentially to zero 
with distance. The Heisenberg model approach predicts 
that, besides the polarization of the system, ferromag- 
netic correlations in the xy plane should be present. In 
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FIG. 8: (Color online) Histogram of the density of the A 
particles as a function of temperature. For large values of Uo, 
one needs very small temperature for the system to polarize 
as the energy associated to the polarization decreases when 
Uo increases. 




FIG. 9: (Color online) The single particle and the anticor- 
related Green functions for U2 > and p — 1 along the x 
axis in the Mott insulator (a) and superfluid phases (b). In 
the MI phase, individual movements of particles are strongly 
suppressed as shown by the exponential decay of Go and Ga 
whereas some anticorrelated movements of particles remain, 
which is shown by the plateau in G a . In the superfluid phase, 
we find a long range order for the Green function correspond- 
ing to the dominant species 0. 

the xy plane spin-spin correlations are measured with the 
correlation function (S Xir +nS Xjr + Sy tr +nS ytr ) which, in 
terms of the particle creation and annihilation operators, 
maps into the anticorrelated Green function G a , Eq. [7] 
In a system where the density is fixed to one particle 
per site due to interactions but where there are differ- 
ent species, there is always the possibility that particles 
move by exchanging their positions. G a measures the co- 
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FIG. 10: (Color online) Histogram of the density of A parti- 
cles for 172 > as a function of t/Uo at fixed total density, 
p = 2. The histogram of the particles is the image of this 
distribution with respect to the line p — 1. The other param- 
eters are L = 8, /3 = 2L/t, U2/U0 = 0.1. The system polarizes 
for t/Uo — 0.05 which corresponds to the MI-SF transition. 
Hence the p — 2 Mott is unpolarized. 



herence at long range of such exchange moves. Indeed 
we see that such exchanges are present in the p = 1 Mott 
phase because G a stays almost constant as distance is 
increased, although it is smaller than its limiting value 
Pop a- This supports the description of the system, at 
strong coupling, in terms of Hcisenberg spins since indi- 
vidual particle degrees of freedom appear to be irrelevant 
to the excitations present in the system. On the other 
hand, spin excitations (i.e. exchanges of particles) ap- 
pear to be relevant. 

In the superfluid phase, the Green function of the dom- 
inant species (0 in the case shown in Fig. [S] (b)) shows 
long range order, indicating the presence of a long range 
phase coherence. We also observe long range coherence 
for the minority species and for the anticorrelated Green 
function. This is typical of a strongly correlated super- 
fluid where different kinds of phase coherence can be 
observed: Phase coherence of the individual particles, 
but also, at the same time, phase coherence of exchange 
moves of particles. 

In the p = 2 MI case, the situation is simpler since 
the phase is not polarized. With positive U%, the system 
will adopt a state where there are two particles of the 
same type on a given site. However, the conversion term 
will couple this state with the corresponding state with 
two particles of the other kind on the site, thus lower- 
ing the energy. On average, those two states have the 
same probability, and the mean number of a given type 
of particle on a site is one. There is also no particular 
density ordering of the system. This behavior is illus- 
trated in Fig. [TOl where the histogram of p\ as a function 
oit/Uo is shown. We see that the transition between the 
unpolarized phase at low t/Uo and the polarized phase 
happens for t/Uo — 0.05. Comparing this with the phase 
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FIG. 11: (Color online) The single particle and the anticor- 
related Green functions for U2 > and p — 2 along the x 
axis in the Mott (a) and superfluid phases (b). In the MI 
phase, individual movements as well as exchanges of particles 
are suppressed, which is shown by the exponential decay of 
all the Green functions. Phase coherence is established in the 
SF phase (b). 



diagram Fig. [3J we find that the polarization occurs pre- 
cisely at the MI-SF transition. Hence the p = 2 MI phase 
is not polarized but the superfluid phase is. 

Neglecting the hopping term, there is no degeneracy in 
the state adopted by the system in the p = 2 Mott phase. 
Then the fluctuations induced by the hopping term do 
not create any phase coherence and all the Green func- 
tions decay exponentially (Fig. QTJa)). In the superfluid 
phase (Fig. ITTT b)). we observe a behavior similar to the 
p = 1 case where all the Green functions reach a plateau 
at large distance, thus showing the presence of differ- 
ent kinds of phase coherence in the system. Due to the 
larger density and the importance of conversion processes 
in this case, the polarization of the system is weaker than 
in the p — 1 case and the leading phenomenon seems to 
be exchanges as G a is larger than Go or G\. 



IV. NEGATIVE U 2 CASE 

In the negative U2 case, the phase diagram predicted 
by MFT, and observed in one-dimensional QMC simu- 
lations is less rich than its positive counterpart. As dis- 
cussed in Section II, we focus on the range of interactions 
\U2\IU0 < 1. Figure [12] shows the MFT phase diagram 
which exhibits MI and SF phases, both of which are un- 
polarized. As U2/U0 — > —1 and for t/Uo = 0, all the 
MI lobes, with even or odd densities, shrink and totally 
disappear at U2/U0 = —1 as the MI phase of density p is 
obtained for (1 + U 2 /U )(p - 1) < p/U < (1 + U 2 /U )p. 
MFT predicts that all phases are unpolarized and all 
transitions continuous. This was confirmed in one di- 



SF unpolarized 




FIG. 12: (Color online) Ground state phase diagram of the 
U2 < case in MFT for p < 4. The surface traces the tran- 
sition points between MI phases (below the surface) and SF 
regions (above). As [C/iz |/C/o — ► 1, the MI regions shrink and 
eventually disappear. 

mension with QMC simulations. Here, we focus on the 
two dimensional case. 



A. Phase diagram and transitions 

Fig. [13] shows the QMC phase diagram for the nega- 
tive U2 case at U2/U0 = —0.1. The boundaries of the MI 
lobes are obtained in the same way as in the positive U? 
case. The results are qualitatively similar to those found 
in one dimension. The system exhibits MI lobes at com- 
mensurate fillings and sufficiently small t/Uo which turn 
superfluid as this parameter increases. In addition the 
system is SF for all incommensurate fillings. We find 
that both these phases are always unpolarized (see be- 
low). As for U2 > 0, the agreement between QMC and 
MFT is better in two dimensions than in one, especially 
at small t/Uo. The agreement is poor as the tips of the 
lobes are approached due to increased quantum fluctua- 
tions. 

In contrast with the U2 > case, all quantum phase 
transitions appear to be continuous in this case. As can 
be observed in Fig. [TJJ there are no signs of possible 
discontinuities in the superfluid density at the transition 
between the MI and SF phases. Other quantities, such 
as p(k = 0) (not shown here) confirm this conclusion. 

B. Polarization and nature of the Mott phases 

To analyse the situation in the p = 1 MI, we again use 
the perturbative mapping to the anisotropic Heiscnberg 
model (Eq. [5]). We restrict ourselves to the case where 
all the interactions are repulsive : — 1 < U 2 /Uq < 0. 
In this range \J xy \ > \Jz\ and the system develops spin 
ordering in the xy plane. In terms of the original model 
representing two types of particles, this spin order in the 
xy plane corresponds to a phase with equal number of 
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FIG. 13: (Color online) QMC phase diagram at low temper- 
ature for three system sizes in the U% < case. Unlike the 
U2 > case, the system never polarizes. The dotted lines 
show the MFT result. 




t/Uo 



FIG. 15: (Color online) Histogram of the density for particles 
of type A for U2 > for a total density p — 1 and U2/U0 = 
—0.1, L — 10 and /3 = 2L/t. The histogram is always centered 
around po = Pa — 0.5 and the system never becomes polarized 
in this case. There is no change of behavior when the MI-SF 
transition is crossed for t/Uo — 0.05. 
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FIG. 14: (Color online) The SF density, p s , as a function of 
t/Uo at fixed density for U2 < 0. The superfluid density varies 
continuously for p = 1 (a) and p = 2 (b), showing no sign of 
a first order phase transition. 



and A particles on each site. This is confirmed by 
QMC simulations. Plotting the density histograms for 
U2/U0 = —0.1, Fig. HH it is seen that the distribution is 
always centered around po = Pa = 0.5 confirming that 
the system is not polarized. In addition, the density- 
density correlations (not shown) do not exhibit any sign 
of density order. 

In this case, the prediction of the mapping onto the 
Heisenberg model is that the system should develop fer- 
romagnetic order in the xy plane and the spin-spin corre- 
lations (S XtT +jiS x , r + S^r+RiS^r) should be long ranged. 
Figure [TO] shows the single particle and the anticorrelatcd 
Green functions in the MI (a) and the SF (b) phases. In 



the p = 1 MI phase, Fig. ITBT a). we see that Go(R) and 
Ga(R) decay exponentially as expected. However, we 
also see that G a (R) quickly saturates to a constant value 
at large separations indicating that exchange moves are 
common. In a system where the number of particles of 
each species is fixed, i.e. no conversion between species, 
this behavior of G a (R) would correspond to a counter 
superfluid phase (CSF) [17| . However, in this system, 
due to the residual effects of the conversion term, it is 
not possible to calculate a corresponding counter super- 
fluid density [ToT ]. However, we do observe the presence 
of a long range coherence of exchange moves as predicted 
by the perturbative mapping to the Heisenberg spin sys- 
tem. The plateau observed in G a takes on its maximum 
possible value at large distances G a — > pop A, showing 
perfect phase coherence, contrary to what was observed 
for U 2 > 0. 

In the superfluid phase at p = 1, Fig. [TC] (b), Go and 
Ga exhibit long range order, indicating the presence of 
phase coherence. The two species remain correlated, of 
course, since the system is still in a strongly interact- 
ing regime (t/Uo = 0.1). This correlation can be ob- 
served, for example, from the fact that G a is larger than 
the product GqG\. This means that, while particles can 
move independently, exchanges of different particles are 
still present. 

The p = 2 Mott phase is once again unpolarized (the 
histogram is similar to the p = 1 histogram shown in 
Fig. n~5|) . In the t/Uo limit, there exist two possible 
degenerate ground states on each site. The first state is 
obtained by putting one particle of each type on a given 
site, thus obtaining an interaction energy Uo + U2 < Uq 
as U2 < 0. The second state is the superposition of a 
state with two particles and a state with two A parti- 
cles: (|00) + |AA))/ v / 2. While the interaction terms Eq. H 
gives in this case an energy Uo, it is reduced to U0 + U2 by 
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FIG. 16: (Color online) The single particle and anticorrelated 
Green functions for JJ% > and p — 1 along the x axis in the 
MI (a) and SF phases (b). For the MI phase, as predicted by 
the spin approximation, we find a dominant G a showing that 
the displacement of particles is due mainly to exchange of 
particles of different types whereas movements of individual 
particles are suppressed. For the superfluid phase, there is 
long range coherence of all the Green functions. 




FIG. 17: (Color online) The single particle and anticorrelated 
Green functions for JJ% < and p — 2 along the x axis in the 
MI (a) and SF phases (b). For the MI phase, there are ex- 
changes of particles due to the degeneracy of the Mott phase 
in the t — limit, whereas individual movements are prohib- 
ited. In the SF, all Green functions show phase coherence of 
individual movements of particles or of anticorrelated move- 
ments. 



V. CONCLUSION 



the conversion term Eq. [31 For these two possible states, 
the mean densities po and Pa are equal and the system 
does not polarize which is directly observed in density 
histograms (not shown here). The p = 2 Mott ground 
state then shows a large degeneracy in the t = on-site 
limit. All these degenerate states are once again coupled 
by second order contributions from the hopping term and 
the degeneracy is lifted by establishing a phase coherence 
of exchange movements. This can be understood in the 
following way: whether a site is occupied by two parti- 
cles, two A particles or one and one A particle, the on 
site energy is the same, because the interaction energy 
is lower for different particles or because it is lowered 
by conversion term for identical particles. Then starting 
from any configuration, exchanges of particles of different 
types will exchange the states of neighboring sites with- 
out changing the on site energy. A phase coherence is 
then established by exchange moves in the Mott phase, 
reminiscent of the one observed in the p = 1 case. This is 
shown in Fig. HPTTa) where we show perfect phase coher- 
ence of the G a function which reaches its limiting value 
of popA = 1 in the Mott phase whereas individual Green 
functions decay exponentially. 

In the superfluid phase, we observe phase coherence 
both for exchange moves and for individual movements 
of particles but, once again, the exchange moves appear 
to give the leading contribution (Fig. UTT b)). 



In this work, we used QMC simulations to deter- 
mine the phase diagram of the spin- 1/2 bosonic Hub- 
bard model on a two-dimensional square lattice. For 
l<2 > 0, we found that the numerical results to be in 
good agreement with MFT, especially at very small t/Uo 
where quantum fluctuations are highly suppressed. The 
system has three phases: A polarized superfluid, an un- 
polarized Mott phase for p = 2 (and likely for all even MI 
lobes) and a polarized Mott phase for p = 1 (and likely 
for all odd MI lobes). We have shown that the polariza- 
tion of this p = 1 Mott phase can be understood in terms 
of an effective anisotropic Hciscnbcrg model and that, 
in the ground state, the phase is completely polarized, 
although extremely low temperatures are needed to ob- 
serve the polarization in the strongly interacting regime. 
The p = 2 phase is not polarized due to the action of 
the conversion term that will transform pairs of identi- 
cal particles into pairs of the other types of particles and 
then suppress possible polarization of the system. The 
first order MI-SF transition for the p = 2 MI phase which 
was predicted by MFT is confirmed numerically. In ad- 
dition, the QMC results showed that in the p = 1 MI 
lobe, moves which exchange the positions of a and a 
A particles are present as evidenced by the saturation of 
the anticorrelated Green function G a . This property was 
not previously addressed by MFT calculations. 

For U2 < 0, the system is always unpolarized and MI- 
SF transitions are all continuous. However, in the MI 
phases, the origin of this absence of polarization is differ- 
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ent for the p = 1 and the p = 2 phases. For p = 1. 
the anisotropic Heisenberg model approach shows the 
domination of an effective in-plane coupling, whereas for 
p = 2, the two possible degenerate states on each site 
are both unpolarized. Here too, we found that particle 
exchange moves are present in the p = 1 and, more sur- 
prisingly, in the p = 2 Mott insulator. 
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